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Abstract 

We evaluate the spin polarization (Edelstein or inverse spin galvanic effect) and 
the spin Hall current induced by an applied electric field by including the weak 
localization corrections for a two-dimensional electron gas. We show that the 
weak localization effects yield logarithmic corrections to both the spin polar¬ 
ization conductivity relating the spin polarization and the electric held and to 
the spin Hall angle relating the spin and charge currents. The renormalization 
of both the spin polarization conductivity and the spin Hall angle combine to 
produce a zero correction to the total spin Hall conductivity as required by an 
exact identity. Suggestions for the experimental observation of the effect are 
given. 
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1. Introduction 

Weak localization (WL) is the result of quantum interference corrections to 
the semiclassical theory of transport [TJ [2]. It manifests itself in good conductors 
as a negative or positive correction to the electrical conductivity depending 
on the symmetry properties of the system. The functional form varies with the 
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effective dimensionality of the sample, behaving as a square root of temperature 
in three dimensions and logarithmically in two dimensions [3]. In the latter 
case, the resummation of the logarithmic correction via the renormalization 
group leads eventually to the Anderson localization transition in d = 2 + e 
dimensions in [S]. In the presence of spin-orbit coupling (SOC), the correction 
is positive and hence manifests as an antilocalizing behavior [6]. SOC affects 
WL because it yields a finite spin relaxation time, which introduces a cutoff in 
the logarithmic singularity associated with the so-called triplet channel of the 
particle-particle ladder, known as the Cooperon. Since the singlet and the triplet 
channels contribute to WL with opposite signs, the elimination of the triplet 
leaves the singlet alone, which then produces the antilocalizing behavior. In 
metallic conductors and doped semiconductors SOC was traditionally attributed 
to the electric field of impurities, which do not affect the nature of the electron 
eigenstates. In the last two decades, however, the two-dimensional electron 
gas (2DEG) has become one of the most analyzed model systems for electrical 
transport, due to the possibility of realizing it in semiconducting systems, and 
more recently at metallic[7] and oxides[8] interfaces. The realization of the 
2DEG leads to the breaking of inversion symmetry with respect to the axis, 
say the z axis, perpendicular to the 2DEG plane, say the x and y plane. In 
these circumstances, in the presence of a finite spin-orbit interaction, Bychkov 
and Rashba have proposed a model Hamiltonian[3], which, besides the standard 
effective-mass kinetic energy term, contains a spin-orbit coupling term linear in 
momentum 

H = ^+a{a'^Py- ( 1 ) 

where p = {Px,Py) is the vector of the components of the momentum operator, 
m is the effective mass and a a SOC constant with and a'^ the standard Pauli 
matrices. The Rashba Hamiltonian Eq.Q has been extensively studied over 
the last twenty years, especially aiming at the development of new spintronic 
functionalities [TU] . In this respect the spin Hall effect fSHEl [TTl IT^ [T?l flTl [T5] 
and the current-induced spin polarization effect [HI [17] (known also as the Edel- 
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stein or inverse spin-galvanic effect) have been the focus of an intensive dedicated 
research. These effects, whose precise definition will be given later on, manifest 
due to the coupling of charge and spin degrees of freedom and hence introduce, 
besides the standard electrical conductivity, new transport parameters. These 
are defined as the linear coefficients relating the spin polarization and the spin 
current to the applied electric field 

((s*/)) = (2) 

where the double brackets indicate the quantum and statistical average, 
and are referred to as the spin polarization or Edelstein and the spin Hall 

conductivities, respectively. As for the electrical conductivity, these transport 
5 parameters can be studied with the well-known impurity technique. One advan¬ 
tage of this technique 1 18]. based on standard diagrammatic perturbation theory, 
is the appearance of the semiclassical Drude-Boltzmann theory of transport at 
the leading approximation in an expansion of the small parameter h/^epr), 
where ep and r are the Fermi energy and the elastic scattering time, the only 
10 two parameters characterizing a disordered Fermi gas. In such an expansion, 
WL arises in the next-to-leading approximation in the expansion in Ti/iepr). 

WL effects in the presence of the Rashba SOC described by Fq.Q have 
been analyzed by several authors, most of the attention having been focused 
on the electrical conductivity only[Tni HOI IHl HH 1131 HU- It is the aim of the 
15 present work to extend this analysis to the other transport parameters men¬ 
tioned above, whose experimental study has developed considerably in the last 
few years [25l [7]. We find that and the spin Hall angle "fsH = drift 
acquire logarithmic corrections which can be absorbed in terms of the renormal¬ 
ization of the scattering time appearing in the electrical conductivity ag. We 
20 emphasize that is not the full spin conductivity which would be 

measured in an experiment |26|. As will be shown in the next Section, can 

be expressed in terms of and The renormalizations of both 

and 'ysH compensate in such a way that has no correction as expected 

on general arguments [27]. 
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25 The plan of the paper is as follows. In the next Section we introduce the 
disordered Rashba model and review the theory of and < 7 ^^^ to the leading 
order in the parameter h/{epT) within the impurity technique. This is necessary 
to prepare the ground for the following Sections. Section III deals with the WL 
localization corrections in the presence of the Rashba SOC. The evaluation of 
30 the electrical conductivity is reviewed as an example. Section IV presents the 
calculation of the WL corrections to and . Section V provides a 

discussion of the results obtained, whereas technical points of the calculations 
are given in the appendices at the end of the paper. From now on, if not 
otherwise specified, we will work in natural units h = c = 1. 

35 2. The disordered Rashba two-dimensional electron gas at leading 

order in l/(eirx) 

In the presence of scattering from impurities, the 2DEG Hamiltonian of 
Eq. 0 acquires an additional random potential term U{r) defined by the aver¬ 
ages 

([/(r))=0, {U{r)U{r')) = (3) 

where r = (x,y) and r' = {x',y') are the coordinate operators, Nq = m/27r 
the two-dimensional density of states and r the elastic scattering time. At 
leading order in the expansion parameter l/{eFT), the selfenergy is given by the 
selfconsistent Born approximation 

where denotes the retarded and advanced Green functions. As discussed 
in [28l [29] , in the presence of Rashba SOG the Green function has a nontrivial 
structure in spin space, whereas the selfenergy remains diagonal, E = 
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G = G°a° + GVi + GV2. Explicitly we have: 

1 


G° 

G^ 

G^ 

G± 


(G+ + G_) 
(G+ - G_) 


2 

Py 

2 


-|^(G+-G_) 


= (e +/X — p^/2to =F ap — E°) 


with 




2t' 


( 5 ) 


( 6 ) 


The Edelstein (EC) and spin Hall (SHC) conductivities are defined in terms of 
the spin polarization and spin Hall current induced by an applied electric field 
taken along the x axis for definiteness’s sake Ej; = —dtA^. The corresponding 
Kubo formulae are 

( 7 ) 


a;—>-0 


and 


CJ 


U) 


( 8 ) 


where the bare vertices = cr*'/2, jy = G^pyl2m and jx = —evx, ^^(p) = 
PxIiTT- — cta^ denote the operators for spin polarization, spin current and charge 
current, respectively. The evaluation of the response functions (7]|8) involves 
the standard bubble diagrams of the Green function lines obtained by the self- 
45 consistent Born approximation Q decorated by the insertion of the impurity 
ladder. This corresponds to the inclusion of the so-called vertex corrections, 
which lead to renormalized vertices[28). 

The expression Q for the EC becomes 


EC 


= -^yTr 

27r 


syG^vx{v)Gi 


( 9 ) 


where the vertex renormalization can be attributed either to the left spin vertex 
or to the right current vertex and we have dropped the dependence on the 
frequency argument of the Green function. In the former case, by using the 
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renormalized spin vertex indicated by a capital letter = ((1 + a;^)/x^)cr^ = 
{tdp = {tdp/T) sy , one obtains [16] 

= —eaNoT, ( 10 ) 


where the subscript 0 in aQ^ indicates the lowest order in the parameter 
l/^epr). We have defined the parameter x = 2apFT and introduced the 
D’yakonov-Perel relaxation time Tpp = 2r((l + pp being the Fermi 

momentum in the absence of the Rashba SOC. The model has two small pa¬ 
rameters l/{eFT) and ajvp with vf = Pf/'^^, in terms of which the above 
paramater x = 4:{a/vF){iFT) can be expressed. 

Similarly, for the SHC one has the expression 


SHC 




( 11 ) 


By using the renormalized spin current vertex + (uj’/2a;)cr^ = + 


55 (vt?/2)(x/( 1 + x'^))Sy [50]. one gets 


-f 


Vf 


2 1 + x^ 


Tr 


syG^Mp)G^ 


^SHC , 


'drift 


+ 


X 


X 


2 l + x^^° 

Vf X 


EC 


EC 


Stt 1 -I- a;^ 2 1 -I- a;^ ° 

^SHC 


( 12 ) 


where we indicated by cF^^ift Hi® SHC in the absence of vertex corrections as 
first computed in Ref. m- It is then not difficult to see that the insertion of 


the result (101 into (121 gives a vanishing SHC. This latter result is actually 

The commutation 


expected following the argument derived by Dimitrova 
relation allows us to write 

ds^ 


dt 


= -2majl, 


(13) 


which, in stationary circumstances, implies ((jf)) = 0. 


The expression (12) shows that the vertex corrections for the SHC are asso¬ 
ciated to the EC. This connection between the two effects acquires a more trans¬ 
parent meaning by adopting the SU(2) gauge-field description of the SOC[3T]. 


6 















In such a picture, the SOC is expressed in terms of a non-Abelian gauge field 
A = A“cr“/2, with Ay = 2ma and = —2ma [32l [33]. The first conse¬ 
quence of resorting to this language is the appearance of an SU(2) magnetic 
field = —(2ma)^, which arises from the non-commuting components of the 
Bychkov-Rashba vector potential. As in the normal Hall effect we have a spin 
Hall drift component of the spin current which can be described as (assuming 

X < 1) 

{{{jy\drift)) = (14) 

with crf^^ = (e/87r)a:^. In addition to the drift current, there is also a “diffusion 
current” due to spin precession around the effective spin-orbit field. Within 
the SU(2) formalism this current arises from the replacement of the ordinary 
derivative with the SU(2) covariant derivative in the expression for the diffusion 
current. The SU(2) covariant derivative, due to the gauge field, is 


VjO = djO + i [Aj,0 ], 


(15) 


with O a given quantity being acted upon. The normal derivative, dj, along 
a given axis j is shifted by the commutator with the gauge field component 
along that same axis. As a result of the replacement 9 —>■ V diffusion-like 
terms, normally proportional to spin density gradients, arise even in uniform 
conditions and the diffusion contribution to the spin current turns out to be 

/7T1 rv 

{{[fyUff)) = (16) 


where D = VpTj2 is the diffusion coefficient. In the diffusive regime the full 
spin current jy can thus be expressed as 


iUy)) 


2ma 


D{{sy)) 


ISH 


(Ux)), 


(17) 


with the spin Hall angle given by 'jsh = and {{jx)) = (^Ex- 

Relation Eq.([T7|, which implies the diffusive limit, coincides with the expression 


(12) when x ^ 1, i.e. when the spin-orbit induced spin splitting is much smaller 


60 than the disorder broadening of the levels. 
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Figure 1: Crossed ladder, also known as Cooperon. The indices a and 5 are 
the incoming spin indices, while /3 and 7 are the outgoing ones. Furthermore, 
the momentum Q is the sum of the initial and final external momenta in each 
Green function line. 


3. The weak-localization correction to the electrical conductivity in 
the presence of Rashba spin-orbit coupling 

In this Section we consider the WL corrections, which arise to first order in 
l/i^epT). In next Subsection we review the evaluation of the Cooperon prop- 
65 agator, whereas in the following one we apply it to the case of the electrical 
conductivity. The spin polarization and spin current response are considered in 
the next Section. 

3.1. The structure of the Cooperon in spin space 

In the absence of SOC, the Cooperon obeys the diffusion equation, which in 
momentum space reads 

where Q = p -|- p' is the sum of the initial and final momenta in each Green 
function line, as shown in Fig. We remind that the particle-hole series 
of the maximally crossed diagrams can be transformed in the particle-particle 
ladder by reversing one of the two Green function lines and the associated 
momentum. The momentum Q corresponds then to the total momentum of 
the particle-particle pair of the Cooperon propagator. The SU(2) gauge-field 
point of view suggests that one can obtain the equation for the Cooperon by a 
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minimal substitution procedure on the momenta of both Green function lines, 
i.e. p — !• p + ^ and p' ^ p' + A , where A is the SU(2) spin-dependent vector 
potentials introduced in the previous Section. The equation for the Cooperon 
becomes then 

2 1 

DjyQ^+AfS‘^rC{Q) = ^^^, (19) 

where 5“ = ct“/2 0 crf/2 is the total spin of the two particles (cr“ 2 refers to the 
upper/lower Green function line). The Cooperon acquires a matrix structure 
acting in the four-dimensional Hilbert space resulting from the combination of 
the two spin one-half spaces. It is convenient to use the triplet-singlet basis 
{I 1t),(l U) + I it))/v^, I U),(l U) - I it))/v^}, in terms of which the 


Cooperon becomes block-diagonal. Since for the singlet the total spin of the 
pair is zero 5 = 0, the one-dimensional corresponding block has the same form 
as in the absence of SOC. The 5 = 1 triplet three-dimensional block can be 


obtained from Eq.(19l by expanding and inverting the square in the left hand 
side. One sees that, besides the standard diffusive terms going as Q^, which 
are diagonal in spin space, linear-in-Q terms arise when considering the double 
product. Finally, the terms arising from the square of the spin-dependent vector 
potential give rise to Q-independent terms. The latter are also diagonal in spin 
space and yield a finite spin relaxation rate. The full expression for the Cooperon 
is obtained by inverting (5“^ (Q) [HI [201 EH E21ES] 


C"^(Q) = 2tvNot 


ie ''^xQ\/Dt 


( DtQ^ + t/tdp 

—ie'^x Qx / Dt DtQ^ + 2t/tdp 


ie-'^'^xQVD^ 


\ 


—ie'^x 


Qy/Dr DtQ^ -\-t/tdp 


\ 


000 Drg 2 y 

( 20 ) 

where 9 is the angle between Q and the x-axis. According to the matrix form 


(20), the singlet state decouples from the triplet states and remains gapless, 
whereas the triplet sector acquires a gap proportional to t/t^p. Furthermore, 
the linear-in-Q terms provide the mixing between the different triplet channels. 
These results are due to the D’yakonov-Perel spin relaxation mechanism, which 
mixes the triplet states and kills their singular contribution. 


We notice that the previous expression (20) is valid in the regime a; ^ 1. In 
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Appendix A we derive the exact form of the Cooperon equation (A.13), whose 


solution tends to expression (20l as a; <C 1. However, the approximate solution 


(20) possesses all the physical features due to the presence of Rashba SOC and 


simplifies the calculations of the quantum corrections. Because of that, in the 


80 rest of this work we will use the Cooperon expression (20) to calculate the weak 
localization corrections. 


3.2. The electrical conductivity to order l/(eFT) 

First we note that the resummation of the maximally crossed diagrams for 
the evaluation of the electrical conductivity gives rise to three diagrams, depicted 
85 in Fig. The following discussion is quite general and applies also to the case 
of EC and SHC to be discussed in the next Section. 

The weak localization correction to the electrical conductivity reads 


Sa — Sa(a) + + (5cr(c) 


( 21 ) 


where the subscript refers to each diagram in Fig. As can be understood by 
looking at Fig. |2a[ the first diagram gives the following contribution 




( 22 ) 


Pfii/X 


Q 


with 


C'^‘^(Q) = -Tr(L'^‘^C'(Q)) 


(23) 


where L^'' = a^ ® is the direct product basis with /x = 0,1, 2, 3, and 


<)(Q)=E9TrG^(-e)^G^ 




m 


K 
J 2 


-Tr 


G§-p( 


, Qx Px qA 


m 


Q-p^ 


(24) 


Notice that the renormalized charge current reads Jx = {—e)px/'rn = (—e)(vx + 


acr^)[55j. In order to evaluate Eq.(24) we make three approximations: first, we 
cut off the sum over Q; second, we expand G(Q) for small Q; third, we neglect 
the weak Q dependence of S'^^^(Q). The natural upper cut-off of the sum over 
Q is provided by the inverse of the mean free path, i.e. the length where the 
diffusive behaviour sets in. On the other hand, the lower cut-off of the sum is 
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Figure 2: Maximally crossed diagrams giving the weak localization corrections 
to the transport coefficients. In each diagram the ladder summation takes into 
account the renormalization of both vertices of the bubble. The crossed lines 
indicate the crossed ladder, also known as Cooperon. In diagram (a) the in¬ 
tegration over the momentum p of four Green functions yields the so-called 
Hikami box, whose full expression is obtained by adding the integrations over 
p and p' of diagrams (b) and (c). 
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given by the inverse of the maximal size of a loop allowed to contribute to the 
coherent backscattering process. Then, for a system of linear size L the lower 
cut-off is 1/L. For small Q the Cooperon is strongly depending on Q and gives 
a divergent contribution which is not modified by the weak Q dependence of 
the weight factors S'^^^(Q). Eventually, we arrive at the following expression for 
the conductivity correction given by the diagram in Fig. 

S , (25) 

p^i/X 

with the weight factors 


Sf. = - V -Tr G^{-e) — a 

(a) 0 p ^ ’ m. 




-Tr 


G_ {-€) — G_ a 


(26) 


and the Q-integrated Cooperon matrix elements 

dO QdQ 


= ^C'^^(Q) = 
Q 


'0 


7i/l 27r 


C^"(Q). 


(27) 


With the approximations introduced above, the contribution (Jcri-j), given by the 
diagram of Fig. is 




(28) 


ppi/X 


with the weight factors 


= -- 


1 


1 


2TrNoT ^ 2 
P,P' 


Finally for 5a(^c) we find 




-Tr 






(29) 


(30) 


ppi/X 


where the weight factors can be expressed in terms of those of 




Tc) 


(31) 


Upon integration over the directions of p, taking into account the Green func¬ 
tion structure of Eq.|5L one discovers that several components of 5”^^^ given 
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in Eq. 
of Eq. 


29 


26) vanish. A similar analysis can be done for the weight factors 

when considering integration over the directions of p and p'. As a 


result, the only non vanishing weight factors are 


cOO _ , 

*(a) - 


\2TrNoT)2DT(^l 


TOP'' 


= e\2T:N^T)2DT 
= e^{2TTNoT)2DT 


4tdp 

3r 

4tdp 


(32) 


and 


00 _ ‘^'pNqT 
4 

2 2ttNqT 




Dt 


X 

(1 + x‘^Y 
^2 


(33) 


Re5(f) = e 


Dt- 


4 ' (1 + x'^Y ’ 

where we have neglected terms beyond the leading order in l/(epT) and the first 


order in a/up. By inserting the weight factors (32) in Eq. (25), the quantum 
correction, given by the diagram of Figj^ is 

o-Q 


<5o'(a) = 


27repT 
2 - 


, L l + a;V2, 

T - 2(iT^0 " 


/2 


4(1 + x"^) 


log 



(34) 


^ / Top ' tIp , j 


The quantum correction, given by the diagrams of Figs. |2b| and [2^ is 


5a, 


(b) 


5a(c) = 



(35) 


To our knowledge the expression for diagrams of Figs. |2b| and [2c] has not been 
given before. In agreement with 


the contribution of the two diagrams 2 b 


and 2c vanishes in the absence of Rashba SOC and does not exhibit scaling 
behavior for a; 7 ^ 0. So it can be ignored and the quantum correction to the 
electrical conductivity is given by 5ai^a) [B ini m m] • 


The above correction (34) can be interpreted as a sum of a localizing contri¬ 
bution from the triplet sector and an antilocalizing contribution from the singlet 
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6(j /(rao) 



Figure 3: The quantum correction to the electrical conductivity as a func¬ 
tion of the variable x = 2apFT for four different values of the ratio Ljl = 
50,250,500,750 from bottom to top. In the y axis we set 5a/[rao), where 
r = l/(27rei?r). 


channel. Therefore, for //L <C a;, the triplet modes acquire a gap that reduces 
their contribution and the total correction changes its sign and becomes antilo¬ 
calizing at a certain x*. The critical value x *, which defines the crossover point 
between WL and WAL, depends on the ratio L/l. In particular, as the ratio Ljl 
gets larger, the value x* decreases pTl IM] . This behavior is evident by looking 
at Fig. for x = 0 the quantum corrections are localizing, whereas for x > x* 
the quantum correction tends to 


Sa = 


ctq , L 


(36) 


which implies the enhancement of the electrical conductivity and the diffusion 
coefficient. 


4. The weak-localization corrections to the EC and SHC 


The equation (21) for the quantum correction to the electrical conductivity 
can easily be generalized to the spin transport coefficients by substituting in the 
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weight factors (26) and (291 the left charge current vertex with the renormalized 


spin density vertex or the renormalized spin current vertex. The correction to 
the Edelstein conductivity then reads 

= + 2RefT';j)C""‘^, (37) 

pfivX 

with weight factors 


=-y-Tr 

(a) o 


^A '^DP ^vr<R^p 


-Tr 

J 2 


^ m ^ 


(38) 


and 


< = -■ 


1 

2'kNqt 


y -Trb^— 

^9 P 2 t p p 


p.p' 


J 2 


-Tr 


G_ G_ ,(-e) —G_ ,cr 
p p m P 


(39) 


Notice that in the right vertex it appears the renormalized charge current vertex 
Jx , whereas in the leading order expression © the bare charge vertex jx appears. 
The non-vanishing weight factors are: 

, -f 2a:® 


r2'|3 


—00 TOP 2ttNqt 2 ' 

IT = e—- :—ar -— 

2r 4 (1 , 

Tjyll Tdp2ttNqT 2 “2;^ + + 2^^ 

VK = —e—- :—ar 

2t 4 

^22 Top 2t:Nqt 

W = —e—- ^—ar 


2t 


(1 -I- x^)® 
2 X^- 1 - 3x'^ -I- X® 

2Y3 


(40) 


(1 -I- x^) 

—33 tdp 2 ttNqt 22;^(1 + 2 x ^) 

W = e— - ^—ar — - 

(l-fx2)3 


2 r 


and = 0 if p A. In the above equation we have defined + 

2ReW^i^. Since the weight factors are diagonal in the indices A and p, only the 
diagonal elements contribute to the correction to the Edelstein conductivity 


Eq.(37). If we take into account only the singlet channel, which is the only mode 
with scaling behavior for x 7 ^ 0 , then all the diagonal elements are related: 


G^^ 


00 


with i = 1 , 2 ,3 and 


GOO ^ 


1 


I 


4t 2 27reFr I ' 


(41) 
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Then, by performing the trace over the product of Pauli matrices in Eq.(37), 
one obtains 


ZTT 


(42) 


from which the quantum correction to the Edelstein conductivity reads 

^ _yy-iQg ^ _ (43) 


^0 1 ^ 

^I ■ 


The above contribution increases the leading order value 

1 . T' 


a^^ = - 




(44) 


2'KepT 

which shows that the logarithmic correction to the Edelstein conductivity can 
be absorbed in terms of the renormalization of the elastic scattering time. Thus, 
Rashba SOC provides the spin relaxation mechanism, which causes the change 
of sign of the weak localization corrections and does not introduce a new scaling 
100 parameter. 

Eor the evaluation of the spin Hall conductivity, we rely on the decomposi¬ 


tion shown in the first line of Eq.(12), where here we replace the bare charge 
current verices jx with the renormalized one Jx and insert a series of maximally 
crossed impurity lines between the two vertices. Hence, the resulting expression 
will have a Cooperon connected to a Hikami box with modified vertices. The 
resulting expression reads 


1 


where 


= -aoSjsH + ■ 

e 2 1 -I- a:"' 


laod^SH = ^ E 


(45) 


(46) 


p^i/X 


Notice that the response function indicated by (ToSjsh contains the renormalized 
charge current vertex Jx in contrast with which has the bare charge 

current vertex. This is the reason why we only find the term (Jq5^sh and we 


don’t find SaojsH in Eq.(45|. The response function aoSjsH describes the 
105 influence of the WL corrections on the spin Hall angle. 
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The weight factors are 


<) = - E jTr 


riAPy ^ nR^p 
i ^ m 2 


-Tr 
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,(-e) — G_ a 
m ^ 


(47) 


and 

- p) 

"('>) 2t:Not 


LPt=-- 
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E T. 


.m2 


-Tr 
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G\G^A-ef-^G\,a^ 


(48) 


The only non-vanishing weight factors are: 

r 2'kNqt 4a;^ -|- 


ZOO ^ _g 


Zii = -e 


m 16 (1-1- x'^P 

T 2'kNqt 2x^ — x® 

m 16 (1-1- x'^Y ’ 


-99 T 2t:Nqt2x 
L = e- 


4x‘‘ 


(49) 


Z33 = e 


m 16 (1-1- x^)o 

r 27riVoT x® 


16 (1 -I- x2)0 ’ 

and = 0 if p A. In the above equation we have defined = 2ReL^^p 
since the contribution of the first diagram of the Hikami boxes is equal to zero, 
^o^lSHXa) = 0- Finally, by replacing the weight factors with the in 


Eq.(42), the quantum correction to the drift component reads 
e x^(l -I- x'^Y 


-aoS'^SH = 


1 

Stt (1-|-x^)o 2'KmD ^ I 


2-KeFT I 


SHC 


(50) 


where we have only taken the logarithimic contribution given by the singlet 


channel. It is then evident that the insertion of the results (501 and (431 into 


(45) gives the vanishing of the weak localization correction to the spin Hall 


conductivity 


Sa 


SHC 


SHC , '^F 


_ I —011 I 

drift ' 




rEC\ . 


as required by the relation (13) in the static limit. 


5. Discussion 

In this Section we develop an interpretation of the structure of the weak 
localization results of the previous Section. By considering Eq. 0 connecting 
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spin current and spin density, we assume for the WL corrections the structure 


<5(0p) = 2ma{6D){{sy)) + 2maD{6{{sy))) + + ^(<5(0.))), 

(52) 

where the variation has been applied to both the coefficients connecting the 
averaged values of the observables and to the average of the observables them- 
110 selves. In the limit cc <C 1, by noticing that VFxj2 = 2maD, we identify the 
second term in the right hand side of the above equation with the second term 


in the right hand side of Eq.(51). In the same limit, i.e. keeping terms up to 
order , we see that the first and last term of Eq. (52) must necessarily cancel 


among themselves. Hence, one can identify the third term of Eq.(52) with the 


115 first term of Eq.(51), thus making clear the reason for introducing the response 
function ctqS^sh describing the quantum correction to the spin Hall angle. 

The direct measurement of the EC or its inverse is possible [3 IM] and a 
possible test of the theory presented here would be the study of the EC as a 
function of an applied magnetic held perpendicular to the 2DEG. On the other 
120 hand, the SHC vanishes in the Rashba disordered 2DEG. It is then important 
to extend this analysis to a case where the SHC is hnite, as for instance when 
also extrinsic SOC from impurities is present [35]. This however is beyond the 
scope of the present work. 

Aknowledgements R.R. thanks Valentina Brosco and Lara Benfatto for 
125 discussions. 


Appendix A. The Cooperon in the presence of Rashba spin-orbit 
coupling 

To describe the sum of maximally crossed diagrams, we have introduced the 
so-called Cooperon, corresponding to the diagram depicted in Eig. In the 
130 main text we have derived an approximate solution to the Cooperon equation, 
by using an argument proposed by Kettemann and Wenk |53| . In this appendix 
we want to obtain the exact derivation of the Cooperon. 
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As can be understood by looking at Fig. the equation for the Cooperon 
can be written as 


C-i(Q) = 27rAfo'r(l - ^ 


(A.l) 


where L^'' = cr^ (g) cr'" and 

/ G« \ 

1 




27rA^oT' 


E 


0 


( G^, ie-^^AG^, -ie^^AG^, 0 ) 


Q-p 


(A.2) 

where fj, = 0,+, —,3, Go = (G+ + G-)j2, AG = (G+ — G-)j2 and (j) = 
aictan{py/px). Since we expect the Cooperon to be large in vicinity of backscat- 
tering p = —p', we set Q = (5(cos 0, sin 0) and expand /(Q) in powers of Q up 
to the quadratic order: 


+/ri)( 0 )Q+> 

where the expansion coefficients are given by 


(A.3) 


2ttNqt 

ir. = ^ 


2ttNqt 


dQ 


Q=0 


(A.4) 


G2) aQ2 


Q=o 


Zeroth order: By symmetry due to the integration over the angle of p, 

there are only three integrals different from zero 7°°^, and . In particular 
one has 

r ..2 ^ 


roo _ 1 

2 (l + :r 2 ) 


= 1 - 


tdp 

/■I— = j~+ — _ — _=_ 

(0) (0) 2(l+a:2) rop 


(A.5) 


where we have introduced the D’yakonov-Perel spin relaxation time r^p. In 
the direct product basis {| ft) ) I ti) ) I it) > I it) } the Q = 0 expression of the 
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Cooperon is: 


C(-; = 27r7VoT(l - /(o)) = 


( t/tdp 0 0 0 ^ 

0 t/tdp t/tdp 0 

0 t/tdp t/tdp 0 

y 0 0 0 t/tdp j 


(A.6) 


First order: At linear order in the Q we have only four contributions: 

(A.7) 


1 P 1 P 


AS 




^ri)W = —= —-/*, 

each one is expressed in terms of the integral 


J~ 


1 


2ttNqt ^ m 


Y. 


G. 


P- 


(A.8) 


where in the Q expansion we have taken the leading order in ajvp and IjepT. 
According to the previous approximation we obtain 


J = 2vpT- 


(A.9) 


(1 + x'^y 

In the direct product basis the linear term in Q has the following expression 


/ 


Cy]{9) = -27rAfor/(i) = 2TrNoT 


TVpX 


2 (l + a; 2)2 


0 ze-*® 


aiS 


0 ^ 

0 ze-*^ 


—ze'"' 0 

—ze*® 0 0 ze 

0 -ze*® -ze*® 


-le 
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(A.IO) 


Second order: The second order contribution is given by five integrals J™) 

, 1(2)1 ^(i) leading order in a/vp and IjepT, these coefficients 
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read 


1(2) = -Dr 


1 - 


x^(6 + 3x^ + 


r2\3 


2{l + x\ 

r++(n^ -2*ep_2;^(6 + 3x2 +x4) 

7 ( 2 ) {0} - -e Dt --- l-'( 2 ) ’ 


(A.ll) 


+_ _ _ x^(6 + 3x^ +x-^) _ , 

A2) on < ^2)3 \2) 


2(1 +x2)3 

where D = v\t j^. In the direct product basis is multiplied for 

( 1 — h{x) 0 0 /i(x)e“^*®/2 ^ 

0 1 — h{x) —h{x) 0 

0 —h{x) 1 — h{x) 0 

y h(x)e^*®/2 0 0 1 — h{x) y 


C'(-J(6>) = 2TrTNoDT 


(A.12) 


where h{x) = x^(6 + 3x^ + x^)/2(l + x^)^. 

Thus, in the triplet singlet basis the Cooperon reads 

2h(x) 


C'-i(Q) ^2TrNoT 


VpT 


DtQ^ + - Sl)(—- DtQ^ 

\tdp 2 

(Q X S) • z + h{x)DTQ^Sl 


(1 + x'^Y 


-DTh{x){Q X S)2- 


(A.13) 


where z is the versor perpendicular to the plane of the 2DEG and S is the total 
spin operator for a pair of electrons. The above result agrees with mm- 
However, the most important effect due to the SOC is the appearance of the 
spin relaxation time, which cut off the triplet channels. This permit us to use 


the diffusive limit Cooperon (20), proposed in 


Appendix B. Integrals of Green functions products 

To perform the calculations of the weak localization corrections of the trans¬ 
port coefficients we encounter the following kind of integral 

^p'Gr(p)Gr(p) (B.i) 

p 
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where ± are the band indices, n, m and I three integer numbers. When the 
chemical potential is the largest energy scale, we are able to transform the sum 
over momentum as 


140 where ^ f2m — p. Then the integral over ^ is carried by evaluating the 
residues to the first order in a/vp and to the zero order in l/{tpT). However, 
this requires some care when the integral involves the two Fermi surfaces and 
the parameter x = 2appT appears. Indeed, since x = {a/vp){4:epT), then when 
X is present one must keep also terms of order l/{epT), which eventually yield 
145 terms of order ajvp when multiplied by x. 

Below we report the results: 

^pGf^(p)G±^(p) = AttNqt^pp[i =f , 

. (B.3) 

J2pG+iP)G-{p)G£(P) = “^ttNot^pp ^ (l T —) , 

(1 =F ^ Up/ 

p 

and those that appear in the Hikami boxes corrected by the impurity line 


^pGf^(p)G^(p) = -2f7riVoT^Pio(^l T , 

P 

^pG^(p)G^?(p)G^(p) = -2iTrNoT^pp ^ , 

p 

^pGf "(p)G^(p) = -2ipNoT^pp , 

G^^(p)G^(p) = -2f7riVoT^(l T , 

P 

G^(p)G^(p)G^(p) = -2*7riVor2^^ , 

P 

^ G^"(p)G^(p) = (l2ixy • 


(B.4) 


Finally, we report the integrals that occur in the weak localization correction to 
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the electrical conductivity 

(p)Gf (p) = dTTiVoT^p^l T , 

^ 1 (B.5) 

^p2G«(p)G«(p)G^(p)G^(p) = AttNot^pI^^ . 

p 
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